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^  ABSTRACT 

The  title  problem  is  treated  in  the  limit  of  large  activation  energy.  It 
is  shown  that  the  evolutionary  process  takes  place  in  a  series  of  distinct 
stages,  and  the  spatial  and  temporal  structure  of  each  stage  is  described.  It 
is  found  that  subsequent  to  thermal  runaway,  the  behavior  of  the  system 
resembles  that  of  self-induced  combustion,  except  that  the  thermal  explosion 
is  now  confined  to  a  thin  surface  layer.  _^-y  4. 
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SIGNIFICANCE  AND  EXPLANATION 

Burning  of  an  initially  cold  combustible  material  is  studied.  It  is 
assumed  that  combustion  is  initiated  by  heating  the  surface  of  the  material. 
This  will  result  in  a  chemical  reaction  first  occurring  at  the  surface,  and 
then  spreading  through  the  material  in  the  form  of  a  flame  or  deflagration 
wave.  Asymptotic  methods  are  used  to  study  the  transient  process  leading  to 
the  establishment  of  the  flame.  Attention  is  confined  to  solid  combustibles 
but  it  is  envisaged  that  the  techniques  developed  here  will  aid  in  the 
eventual  understanding  of  the  evolution  of  gaseous  combustion. 


EVOLUTION  OF  DEFLAGRATION  IN  A  COLD  COMBUSTIBLE 


SUBJECTED  TO  A  UNIFORM  ENERGY  FLUX 
A,  K.  Kapila* 

1.  INTRODUCTION 

Consider  a  combustible  material  confined  to  a  container  and  capable  of  underarm,: 
exothermic  chemical  reaction  of  the  Arrhenius  type.  There  are  essentially  two  modes  by 
which  burning  can  be  initiated  in  this  system:  either  throuah  self-induced  heating  (th^vma'. 
explosion)  or  by  an  external  stimulus  (ignition).  In  the  first  mode  the  initial 
temperature  of  the  material  is  high  enough  so  that  a  significant,  amount  of  chemical  Vat 
release  begins  to  occur  immediately  throughout  the  system.  If  heat  loss  is  not 
overpowering  the  temperature  will  rise,  but  seldom  uniformly  in  space.  The  fastest  rise 
usually  occurs  at  one  or  more  discrete  sites  where  hot  spots  form  and  eventually  combust  ion 
waves  originate. 

A  detailed  mathematical  treatment  of  the  dynamics  of  self-induced  combustion  for  a 
dif fusional-thermal  model  has  recently  been  given  by  Kapila  [1].  A  symmetric  geometry  is 
chosen  so  as  to  lead  to  the  development  of  precisely  one  hot  spot,  and  the  spatial  and 
temporal  structures  of  the  various  phases  of  the  evolutionary  process  are  described.  It  is 
shown  that  while  the  birth  of  the  hot  spot  is  a  rather  slew  phenomenon,  its  growth  and 
transition  into  a  deflagrating  wave  are  extremely  rapid  events  indeed.  A  partial 
development  alonq  the  same  lines,  undertaken  independently  and  including  extensions  to 
several  geometries,  is  due  to  Fassoy  and  Poland  [2-4]. 

The  above  treatments  have  all  been  asymptotic,  based  upon  the  limit  of  large 
activation  energy.  This  limit,  although  in  implicit  use  since  the  early  work  of  Zeldovich 
and  Frank-Kamenetski i  (see,  for  example,  [5]  and  [6]),  was  suggested  as  the  basis  for  a 
rational  asymptotic  theory  only  recently  by  Williams  [7).  Since  then  it  has  been  applied 
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successful ly  to  a  wide  variety  of  problems  in  combustion  theory,  and  in  fact  has  inspired  a 
forthcoming  monograph  on  the  subject  by  Buckmaster  and  Ludford  [«].  In  the  present 
context,  Kassoy's  extensive  investigations  (see  [9]  and  the  references  contained  therein) 
of  the  spatially  uniform  thermal  explosion  deserve  snecial  mention. 

In  the  ignition  mode  the  initial  temperature  of  the  combustible  is  so  low  that  it  is 
practically  inert,  and  will  not  burn  if  left  alone.  Enerqy  transfer  from  an  external 
source  is  therefore  needed  to  activate  the  chemistry.  This  raises  a  natural  ouestion:  how 
long  should  the  external  stimulus  be  maintained  before  the  material  is  sufficiently  excited 
to  develop  self-sustained  combusion?  Attempts  to  treat  this  question  embody  the  classical 
theory  of  ignition.  For  solid  combustibles  the  status  of  the  theory,  and  the  pertinent 
literature  as  of  1^70,  are  contined  in  the  review  article  by  Merzhanov  and  Averson  [in]. 
From  a  theoretical  standpoint,  two  problems  appear  to  have  received  the  most  attention. 
These  treat  a  homogeneous,  semi-infinite  reactive  solid,  heated  either  by  applying  a 
constant  energy  flux  at  the  surface  or  by  raising  and  subsequently  maintaining  the  surface 
temperature.  The  nonlinear  parabolic  differential  eouations  governing  the  problems  have 
either  been  solved  numerically  or  by  a  variety  of  approximating  techniques  [10].  They  all 
have  one  common  feature:  ignition  is  deemed  to  have  occurred  when  the  solution  is 
perceived  to  have  departed  significantly  from  the  chemically  inert  solution  for  the  same 
configuration. 

The  large  activation-energy  limit  has  played  a  clarifying  role  in  ignition  theory  as 
well.  Thus,  most  past  approximations  were  rendered  obsolete  by  the  recent  asymptotic 
studies  of  Linan  and  Williams  [11,  12].  Their  vrork  has  two  important  features.  First, 
the  onset  of  ignition  is  defined  by  a  mathematically  precise  criterion  based  upon  "thermal 
runaway".  In  other  words,  the  solution  is  sought  as  a  small  deviation  from  the 
corresponding  inert  solution,  and  the  instant  of  ignition  is  defined  to  be  the  one  at  which 
the  perturbation  develops  a  mathematical  singularity.  Second,  explicit  asymptotic 
expressions  are  given  for  the  ignition  delay  period.  The  results  are  found  to  be  in  close 
agreement  with  numerical  values  available  in  the  literature. 


While  the  above  work  answers  one  rruestion  satisf actor ilv ,  it  raises  another,  nanelv, 
what  is  the  post-ignit ion  transient  that  culminates  in  the  development  of  a  def  laarat  ion 
wave?  A  description  of  the  transient  falls  outside  the  scope  of  ianition  theory,  since 
large  departures  from  inert  heating  will  be  involved.  The  singularity  of  Linan  and 
Williams'  quasi-steady  solution  does  provide  a  clue,  however;  it  sucrctests  that  unsteady 
effects  will  have  to  he  reinstated  for  further  development  of  the  solution. 

This  paper  aims  at  describing  the  entire  seouence  of  events,  subsenuent  to  ignition, 
for  the  case  of  the  constant  surface-flux  stimulus.  We  shall  find  t^at  the  situation  is 
quite  akin  to  that  for  the  self-induced  combustion  mode  discussed  in  r  1 '  ,  except  tha^  th*j 
thermal  explosion  is  now  confined  to  a  thin  surface  layer.  The  analysis  ana in  Iraws  upon 
large  activation-energy  asymptotics,  and  a  solution  uniformly  valid  in  space  an \  time  is 
sought.  The  problem  corresponding  to  a  step  increase  in  surface  temperature  will  be  the 
subject  of  another  paper. 

On  physical  grounds  it  is  nuite  clear  that  for  a  solid  reactant  the  ultimate  state 
will  be  one  of  nearly  steady  deflagration  if  heat  loss  is  absent.  For  a  gaseous  mixture 
the  situation  is  not  so  obvious.  Depending  upon  the  strength  of  the  reaction  and  its 
coupling  with  the  qas  dynamics,  the  final  state  may  well  be  a  detonation.  It  is  expected 
that  the  techniques  developed  here  will  aid  in  the  eventual  under stand  inn  of  t^e  evolution 
of  combustion  in  confined  gases. 


2.  MATHEMATICAL  FORMULATION. 

Consider  a  combustible  material  occupying  the  right-half  space  x'  >  0  and  capable  of 
undergoing  a  first-order,  one-step  Arrhenius-type  chemical  reaction.  (Primes  indicate 
dimensional  quantities.)  We  envisage  a  diffusional-thermal  model  in  which  the  material  is 
assumed  to  have  constant  physical  properties  (density  p',  thermal  conductivity  k' , 
material  diffusivity  D'  for  the  reactant  species,  and  specific  heat  c').  Initially  the 
temperature  is  0J  and  the  fuel  mass  fraction  Y^.  Commencing  at  time  t'  -  0,  a 
constant  heat  flux  q'  is  applied  at  the  surface.  The  equations  governing  the  thermo¬ 
chemical  behavior  of  the  material  are 


p'ce;, 

-  k'9x.x.  +  Q'A'Y' exp( 

-E’/R'e')  , 

p’Yt> 

-  P'D'YX.X.  *  *'Y'exp< 

-EVR'8’)  , 

(0,f) 

-  -q’  ,  e*(«,f  )  -  8J; 

Y’,  (0,f  )  -  0,  Y*  (m,t 1  ] 

1  -  1 

© 

X 

<D 

-  0',  Y'(x' ,0)  -  Y’  . 

Here  Q'.E'.A*  and  r'  are,  respectively,  the  heat  of  reaction,  the  activation  energv, 

the  pre-exponential  factor  and  the  gas  constant.  Cn  selecting  0^,  Y^ ,  k'0^/q'  and 
2  2 

p'c'k'0^  /q'  to  be  the  units  of  temperature,  mass  fraction,  length  and  time 
respectively,  the  above  system  reduces  to  the  following  dimensionless  form: 


0fc  -  6xx  +  AY  exp(-E/0)  , 

Yfc  »  iT’y  -  (A/B)Y  exp( -E/0 )  , 


8x(0,t) 


Yx(0,t)  -  0,  Y(»,t) 


-1,  9<-,t)  -  1  , 

1  , 

0(x,O)  »  1  , 
y(x,0)  -  1  . 

The  nondimensional  parameters  appearing  above  are 

A  ”  (fi'A'Ypk'0^/q’ 2,  the  Damkohler  number, 

B  ■  (9'Y^)/( c'0^ ) ,  the  heat-release  parameter, 

L  -  k'/tp'c'D'),  the  Lewis  number,  and 
E  «  E'/(R'0p),  the  activation  energy. 


(2.1a) 

(2.1b) 

(2.2a) 

(2.2b) 

(2.3a) 

(2.3b) 
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The  mathematical  model  considered  above  is  quite  adequate  for  an  approximate 


description  of  solid  combustion  (but  not  for  gaseous  reactants,  since  material  deformation 
has  been  ignored)  and  becomes  still  more  realistic  by  taking  L  =  “,  because  diffusion  of 
the  solid  reactant  relative  to  that  of  heat  is  negligible.  We  shall,  however,  present  the 
analytical  details  only  for  L  =  1;  the  special  symmetry  of  the  equations  then  renders  the 
mathematics  particularly  simple  without  compromising  the  essential  physics.  This  was 
verified  by  actual  computations  for  L  =  “,  and  it  was  found  that  the  asymptotics  must 
then  be  taken  to  higher  orders  to  capture  the  eauivalent  results.  A  discussion  of 
significant  departures  from  the  L  =  1  case  is  given  in  the  last  section. 

Although,  as  pointed  out  above,  the  model  under  consideration  is  deficient  as  a  tool 
for  studying  the  evolution  of  gaseous  combustion,  it  has  provided  satisfactory  descriptions 
of  several  other  aspects  of  gaseous  flames.  See,  for  example,  accounts  of  cellular  flames 
[13]  and  of  flames  in  accelerating  environments  [14];  also  see  [15]  for  a  rational 
derivation  of  the  model  from  the  equations  of  aerothermochemistry . 

Typical  values  of  A, B  and  E  are  lO1^,  2  and  25,  respectively  (see,  for  example, 
Bradley  [16]).  It  is  particularly  realistic,  therefore,  to  attempt  an  asymptotic  analysis 
in  the  limit  E  +  ■»,  taking  B  to  be  an  0(1)  parameter.  The  largeness  of  A  is 
accounted  for  by  setting 

A  =  A  exp(E/0c)  (2.4) 

where  A  is  assumed  to  be  at  most  algebraic  in  E  for  large  E,  and  will  be  assigned 
during  the  course  of  the  analysis.  Thus  (2.4)  defines  6c,  which  is  assumed  to  be 
0(1).  Also,  we  let 

1  <  8  <  1  +  B  .  (2.5) 

c 

The  lower  limit  ensures  that  combustion  is  not  instantaneous,  while  the  upper  limit,  as  we 
shall  see,  allows  a  self-sustained  deflagration  wave  to  propagate  into  the  cold  material. 

For  unit  Lewis  number,  equations  (2. 1-2.3)  show  that 

8+BY=0I+B,  (2.6) 

where  fl  ,  the  inert  temperature,  is  the  solution  of  (2.1a),  (2.2a)  and  (2.3a)  for  A  = 


0,  and  is  given  by 


1  ♦  2  f  t  /IT  ) 


t  r.  ” 


The  result  (2*6)  effectively 
variable.  Following  Linan 
where 


1/2  2 

^xn'-x  /(4tP  -  x  prfc'x/f  ?►  t- V  • 
eliminates  Y  from  explicit  consi ierat ion  as 
and  Williams  f 10?  we  recast  the  0 -problem  in 


a  dependent 
t  errrs  c  f  *  , 


*  =  0  -  o  f  **  Q  ) 

I 

is  the  departure  of  9  from  the  inert  9^.  Equations  (2.4),  (2.^)  and  (2.P)  then  lead  to 
the  reduced  system 

$  =  4>xx  +  (A/UHR  -  $  )exp?  f  O^/c  ){6’’  -  (9  *  s'1)}’  ,  I7.“l 

♦  (0,t)  =  i(«,M  =  o  ,  12.  ID) 

X 

<t  (x,0)  =  0  ,  f  2  .  11) 

where 

£  =  92/F  .  12. 1.1) 

c 

Our  aim  is  to  develop  an  asymptotic  solution,  uniformly  valid  in  x  and  t,  as  n  ♦  n. 

The  solution  will  be  seen  to  evolve  in  several  staaes,  each  character i zed  by  a  distinct 

time  scale.  The  longest  of  these  is  the  inert  stage,  of  0(1)  extent  on  the  time  scale 

t  and  characterized  hv  essentially  zero  chemical  activity.  A  short  transition  period 

of  0(e)  duration  follows,  in  which  a  weak  chemical  reaction  develops  in  a  tMn  surface 

layer,  eventually  leading  to  thermal  runaway.  The  reaction  intensifies  during  the  ignition 
2 

state,  0(e  )  long  in  t,  and  a  hot  spot  appears  at  the  surface.  The  hot  spot  shrinks 
and  strengthens  in  an  exponent  .ally  brief  explosion  period.  A  still  shorter  detachment 
stage  follows,  in  which  the  hot  spot  leaves  the  surface  and  heains  its  advance  into  the 
interior.  Finally,  a  well-defined  deflagration  wave,  traveling  with  practically  uniform 
speed,  is  established  during  the  slightly  longer  propagation  staoe.  Thus,  in  physical 
terms,  the  transient  beyond  the  transition  staoe  is  similar  to  that  in  self- induced 
combustion,  except  that  the  thermal  explosion  is  now  confined  to  a  narrow  surface  layer. 

We  shall  take  advantage  of  this  similarity  to  omit  details  whenever  possible,  referring  the 
reader  to  [1]  instead.  In  the  interest  of  clarity,  however,  some  repetition  has  been 
unavoidable. 

A  schematic  view  of  the  transient  appears  in  Figure  1. 
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INERT  AND  TRANSITION  STAGES 


The  first  two  stages  were  treated  in  detail  by  Linan  and  Williams  [11].  In  order 
continue  the  solution  further  it  is  helpful  to  briefly  reconstruct,  and  slichtly  amplify, 
their  analysis. 

For  small  t,  0=0  +  $  <  0  ,  so  that  the  reaction  term  in  (2.9)  is  exponential'.*.- 

I  c 

small.  Therefore,  the  solution  to  (2.9-2.11)  is  inert,  i.e. 

4>  =  0  +  e  st  ,  '  3  .  1  ' 

where  est  stands  for  exponentially  small  terms.  Correspondinaly ,  (2.*)  shows  that 

Y  =  1  +  est  , 

i.e.,  reactant  consumption  has  been  negligible.  Itiis  solution  becomes  void  and  aives  way 

to  the  transition  stage  as  soon  as  the  chemical  term  comes  into  play  in  (2.9).  This  will 

occur  first  at  the  surface  (where  temperature  is  at  its  maximum)  at  time  tc  such  that 

6(0,t  )=  0  ,  i.e.  from  (2.7), 

I  c  c 

t  =  tt(0  -  1)2/4  .  (3.2' 

c  c 

The  argument  of  the  exponential  in  (2.9)  indicates  that  for  £  *  0,  chemical  activity  wil 
be  confined  to  the  region  where  9_  +  -  6  «  0(e),  i.e.  to  the  thin  surface  zone  R, 

(see  Figure  2)  where  the  appropriate  independent  variables  are  £  and  x,  defined  by 


3$  /3£  -  0  at  £  =  0;  i  =  0,1,2,... 


while  the  initial  conditions 


$.(£•-»)  =  0 
l 

a  re  a  result  of  matching  with  the  inert  solution  (3.11*  ^e  *irst  two  terms  m  m,*'  a  re 
found  to  be  (see  [11]) 


' 


=  An(f>  +  P  )rfT)  +  f ,  <  t  ) 


where 


P(T)  =  exp r  f  (t  )  +  (*t_) 


- 1  n 


and  the  yet  unknown  functions  satisfv 


f p(-«)  =  f  (  —  )  =•  o 


(  3.Q  ' 


As  Z  +  equations  (3.f<)  show  that  do  not  d^rav  to  zero,  thereby  violatinr 

the  right  boundary  condition  in  (2.10).  This  suggests  the  presence  of  a  new  reaior. 
separatina  R^  from  the  inert  region  Rp ,  as  shown  in  Figure  2.  In  tue  spatial 

coordinate  is  taken  to  be  X,  defined  bv 

x  =  /e  x  or  %  =  X/Ve,  X  >  o  , 

which  reduces  (2.9)  to 

9  =  +  /c  (A  /R)(B  -  4  )exp!  (6?/c  )  {5  '-(6  -  /e  X  *  i  ♦  o  |/e>)  '}•  . 

Anticioating  to  be  0(e),  it  is  clear  that  the  reaction  term  in  the  above  emiation  5  s 

exponentially  small.  Oi  setting 


3/0 

=  ei>0(X,t)  +  e  tyx.T)  +  ...  , 


(  3 . a  1 


ip^  and  ip  ^  are  found  to  satisfy 


0  =  iin  ;  T  i  —0° .  X  >  n  , 

t  XX 


\>n(0,T)  =  fp(T),  4>0  (0,T)  =  -AnP(T),  l|'nf”,Tl  =  h  ,  f 


<3. 10) 


)>0(X,-»)  =  n  , 


and 


►-I  “  <•  T  >  -  X  >  0  , 

t  XX 


(  0  ,T  )  =  f  ^  (  T  )  ,  lp1  (  0  ,  T  )  =  (3/2)A^!P(T)l2  -  A()P(T)f1(T)  ,  (“,T  ) 

^(X,  — )  -  0  . 


•} (?- 


n ) 
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The  left  boundary  cond  it  ions*  in  t-he  two  problems  above  result  fron  mat  "rhino  with  R-,, 


while  the  right  boundary  conditions  and  the  initial  conditions  are  provided  by  "'atchir.o 
with  the  inert  region  Rg.  Both  the  prob.ems  can  be  solved  by  usina  a  Larlace  transform  in 
n.  Thus,  (3.10)  has  the  transformed  sol  tit  ion 


t^0 ( x  ,T )  *  L.T.  of  ^0(X,T  ) 

T  2 

*  /  r-sf  (o)  +  a.  exp{f_(o)  +  (wt  )  1,/2o}:e‘  l' 
J  0  V  0  c 


f  3, 12  ) 


do  . 


An  integral  equation  for  the  still-unknown  function  fg(i)  can  be  obtained  by  invert  mu 
^  (s,T)  and  taking  the  limit  n  0.  It  is  convenient  to  write 

VT)  =  '  (3.13) 

where 

T  «  In [ ( fl t  )V4AJ  +  (lit  )'1/2I  .  '3.14) 

C  0  c 

Then  fQ  satisfies  the  parameter-f ree  integral  equation 


f0(t, 


(TT  ) 


•1/2 


/  <3 


o)  1/'2exp[f  Q(o)  +  a]  do 


Linan  and  Williams  [11]  obtain  an  equivalent  integral  equation,  solve  it  numerically  and 
show  that  fg  is  positive,  monotone  increasing  and  becomes  unbounded  at  i  -  b^ ;  they 
find  to  have  the  numerical  value 

bQ  =  -0.431  . 


Following  [11]  we  choose 


An  =  tut  )_1/4e  °  (3.15) 

u  c 

which,  through  (3.13)  and  (3.14),  places  the  singularity  in  f^fi)  at  t  =  0 .  The 
following  asymptotic  expansion  can  be  developed: 


♦Specification  of  (0,r)  in  (3.11)  requires  carryino  (3.S)  to  an  additional  term. 
X 
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f  0  f  T )  =  -(1/2)ln(-r)  +  ln[/it/(2A0  )]  +  o(  1  >  as  T  ■*  n  .  '-.i 

(The  above  expansion  also  appears  in  [11]  but  there  the  second  term  seems  to  have  a 

misprint.)  Once  fp  is  known,  (3.12)  can  be  inverted  to  yield  i(i  (X,t).  We  shall  nee4 
the  behavior  of  i^fX/CI)  for  small  X,  and  this  can  he  calculated  from  the  large  s 
expansion  of  ^  (s,0).  The  result  is 

^p(X,0)  =  -In  X  +  [ln{/w/( 2Ap ) }  -  y/2]  +  0(X2)  as  X  ♦  0  ,  (3.1 

where  y  is  the  Euler's  constant. 

The  unknown  function  f^(i)  can  also  be  determined  from  (3.11)  in  an  entirely 
analogous  way,  and  is  found  to  have  the  asymptotic  behavior 

fjd)  =  -(3/4)7T3/2(2  -  IT  )  *’ (-T  )*1/2  +  o(1)  as  T  *  0_  . 

Returning  now  to  the  region  Rj,  we  note  that  the  unboundedness  of  fp  and  f, 

leads  to  a  breakdown  of  the  expansion  (3.5)  as  t  ♦  P  .  This  is  as  far  as  Linan  and 
Williams  [11]  went;  they  identified  the  appearance  of  the  singularity  as  thermal  runaway 
i.e.  the  onset  of  significant  departures  from  the  inert  solution. 
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4.  IGNITION  STAGE. 

In  order  to  continue  the  solution  further  in  time,  we  observe  that  (3.5)  yields 

~  1  £  <_T  *  1  as  T  ■*  0  ' 

indicating  that  the  hitherto  neglected  unsteady  term  in  (3.3)  must  be  reinstated  when 
T  =  0(e).  This  suggests  a  new  time-variable  T,  defined  via  the  scaling 

2 

t  =  £T ,  or  equivalently,  t  =  e  T. 

As  Figure  2  indicates,  the  relevant  regions  in  the  xt-plane  are  F^  (characterized  by 
£  and  T)  and  R?  (characterized  by  x  and  T).  Of  these,  the  latter  is  easily 


disposed  of;  in  it  (2.9)  reduces  to 


»T  =  E*XX  *  6St- 


Therefore,  to  leading  order,  $  can  depend  only  on  X.  Matching  with  R1  then  gives 


♦  *  ei(io(x,0)  +  o(e ) 


In  R^,  (2.9)  assumes  the  form 


♦  T  =  +  e3/2(A0/B)(B  -  $  )exp[(92/£ )  {e^1  -  Oc  +  4>  -  e£  +  o(e))  ’}, 

£  >  P,  T  >  -»  . 

Matching  with  R^  as  T  ♦  -  «  dictates  the  expansion 

*  *  ■  lne  ♦  e(01(£,T)  -  In  aq]  +  ...  , 

which,  when  substituted  into  (4.2),  yields 

U  =  U  +  exp(u  -£),  £  >  0,  T  >  -  “  . 

T  1££ 


The  left  boundary  condition  in  (2.10)  provides 


U  (P,T)  =  0 

£ 


while  matching  with  the  R^-expansion  (4.1)  yields  (cf.  (3.17)) 


In  £  +  (In  ^  y )  as  £ 


The  initial  condition 


U.  ♦  -  4  ln(-T)  +  In  ^  -  -T—z  v„(C)  as  T  *  -  “  I4.4-1! 

is  the  result  of  matchinq  with  the  ^-expansion  (3.5)  •  Here, 

V«>  =f  ,2(2  -  ,)'1  +  ^*<5  ♦  .-6)  .  M.- 

We  observe  first  that  the  above  problem  for  is  parameter-free,  and  therefore 

defines  a  universal  function.  Second,  each  of  the  three  terms  present  in  the  full  equation 
(2.9)  still  survives  in  the  reduced  form  (4.4a),  albeit  the  Arrrhenius  exponent  has  been 
simplified  via  linearization  about  a  reference  state  (i.e.  the  inert  state).  This  feature 
is  reminiscent  of  the  induction-period  problem  for  self-induced  combustion,  discussed  in 
[1),  and  we  shall  find  that  the  similarity  extends  to  the  nature  of  the  solution  as  well. 
Equations  (4.4)  were  solved  numerically  by  spacewise  discretization,  followed  by  a 
numerical  integration  of  the  resulting  system  of  ordinary  differential  equations  with  the 
stiff  ODE  solver  EPISODE  B  [17).  Prior  to  discretization,  a  coordinate  transformation  was 
used  to  map  the  problem  onto  a  finite  spatial  domain. 

The  numerical  solution  is  displayed  in  Eioure  3.  Initially  the  temperature  rise  is 
everywhere  gradual,  but  then  t?^  begins  to  increase  rapidly  near  F,  =  0  while  variations 
are  more  leisurely  elsewhere.  Eventually,  at  a  definite  time  T—  (=  -ia.719),  n  (0, T) 
becomes  unbounded.  The  evolution  of  the  singularity  at  T  can  be  examined 
analytically.  (The  following  discussion  will  he  brief,  since  the  situation  is  analogous  to 
that  in  [1).)  We  let 

T  =  Too-p,  p>0,  (4.6) 

and  in  the  limit  p  ♦  0,  develop  two  asymptotic  expansions  for  U^.  An  outer  expansion, 
valid  awav  from  5  =  0,  has  the  form 

U,  ”  U,0(t)  *  P  Ull(5)  +  '  (4.7) 

where  U1()(5)  is  the  numerically  computed  solution  of  (4.4)  at  T  =  T^,  and  the  higher- 
order  terms  0^  can  be  determined  by  substituting  (4.7)  into  (4.4a).  An  inner  region,  or 

boundary  layer,  is  characterized  by  the  variables  p  and  n,  where 

h  =  C/Zp,  0  <  n  <  00  ,  (4.8) 

and  the  inner  expansion  is  seen  to  have  the  form 
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ot  *  -  In  p  +  o  (n)  ♦  /p  q  (n)  ♦  •••  .  4,  ' 

In  the  limit  p  ♦  0,  n  fixed,  the  boundary^layer  thickness  is  0(/t),  and  th*»  expand*.  • 
(4.9)  indicates  explicitly  that  the  singularity  in  is  logarithmic.  When  (4.°)  is 

substituted  into  (4.4a),  we  obtain  the  following  equation  for  gp: 

g£  -  ~  n  g^  +  exp  g  -  1  *  0#  0  <  n  <  00  •  t4.1pa 

The  condition 

g’(0)  =  0 

is  by  now  familiar,  while  the  behavior 

-2 

=  -2  In  n  +  a^  +  o ( n  )  as  n  *  " 

comes  from  suppressing  the  exponential  growth  of  gQ  for  n  ♦  »;  such  a  growth  wouli  no* 
be  commensurate  with  the  slow  rate  of  change  observed  outside  the  boundary  layer.  The 
constant  aQ  is  determined  by  matching  the  expansions  (4.7)  and  (4.9),  which  reouires 

U1fl(C)  =  -2  In  C  +  aQ  +  0(5)  as  5*0.  (4.11' 

A  comparison  of  the  above  asymptotic  expression  for  U  with  the  numerical  solution  of 
(4.4)  at  the  edge  of  the  boundary  layer  (i.e.  for  such  a  choice  (£,T)  which  satisfies 
£  «  1,  T  -  T  «  1  and  n  »  £/Vt  -  T  »  1)  then  yields  a  (*  2.5).  With  a„  known, 

ao  oo  U  U 

gQ  can  be  computed  by  integrating  (4.10)  backwards.  The  result  is  a  monotonic  curve  with 

gn(0)  >0  (4.12) 

(see  Figure  5  of  (1).  The  higher  order  g^  occurring  in  (4.9)  can,  in  principle,  be 
computed  in  an  analogous  way. 

To  summarize,  we  have  shown  that  the  ignition  stage  terminates  at  a  definite  time 
T—.  For  T  <  T  ,  the  tempei  ure  increases  throughout  the  0(e)  thick  surface  layer 

but  remains  an  o(e  In  e)  perturbation  of  the  inert  temperature  (cf.  (4.3)).  Toward 
the  end  of  the  ignition  period  (T  *  T^  or  p  ♦  0),  a  boundary  layer,  or  hot  spot,  of 
thickness  0(/p  e)  develops  within  Rj.  In  the  hot  spot  the  solution  grows  rapidly  and 
has  the  expansion 

♦  “  £  In  e  *  e  (-In  p  +  <gQ  ( n  >  -  In  aq  +  0(/p) )  +  o(e),  n  >  o  ,  (4.13) 

while  in  the  remainder  of  R^  the  solution  remains 
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(4. 14  ) 


♦  *  -  J  e  In  E  +  £[U  (C)  -  In  AQ  +  0(p>)  +  o(E),  C  >  0  .  (4.14) 

In  the  0(/e)  thick  region  R1  the  solution  is  constant  in  time,  to  leading  order,  durina 
the  entire  ignition  stage  (cf.  (4.1))  while  in  the  inert  region  ,  $  =  0  to  all 

algebraic  orders  in  e . 

As  the  process  evolves  further,  0(1)  departures  from  the  inert  solution  will  begin 
to  occur  at  the  surface.  Therefore,  it  will  be  convenient  to  discard  $  in  favor  of  9 
as  the  dependent  temperature  variable. 


FVPLOSI0N  STAG*1. 


5 , 

As  p  ♦  b,  the  0(/p  f.  )  hot  spot  diminishes  in  size  hut  orows  in  intensity, 
eventually  causinq  the  expansion  (4.12)  to  break  down  when  e  In  p  =  0(1).  The  further 
development  of  the  solution  therefore  takes  place  on  the  exponentially  short  tine  seal* 
o,  defined  by 

o  =  -  e  ln(/e  A  p)  ,  o  >  0  and  r>M  )  . 

0 

correspondingly,  £  is  exponentially  snail,  and  in  particular, 

e  in  f,  ^  -  0(1) 

at  the  edge  of  the  hot  spot . 

The  analysis  again  proceeds  as  in  [1]  and  only  a  brief  description  will  he  diver. 

here.  In  terns  of  the  new  hot  spot  variables  n  and  o,  equation  (2.Q)  reads 

^  n6  -  e6  +  (e/BXB  +e  -  6)expt(82/e){e*2(8  -  o)  -  6  ~ 1  > '  (S.rM 

nn  2  n  o  i  ccc 

where  (2.R)  has  been  emploved  to  renlace  f  by  9.  Fron  (2.7)  and  (3.2)  we  note  that  in 
the  hot  spot, 

°1  =  6c  +  est  *  (5.3) 

Also,  the  left  boundary  condition  in(2.2a)  reduces  to 

30/9n  »  o  +  est  at  n  =  n  ,  (^.4) 

i.e.  the  time  scale  of  the  explosion  stage  is  so  small  that  practically  no  heat  enters  the 
surface.  It  can  then  be  shown  that  the  solution  of  (5.3)  satisfying  (5.4)  and  matching 
with  the  ianition  stage  has  the  expansion 

,  .  -  (8  -0 ) 2 { 0  +B-62(6  -o)'1) 

8  »  8  (6  -o  )  +  €0  (0  -o)  [o  (n)  -  In  — - - - - — — - ]  +  ...  .  (8.5) 

c  c  c  c  0  2 

P8 

c 

Here  g^fn)  is  the  same  function  that  appeared  in  (4.°),  i.e.  the  hot  spot  retains  the 

spatial  structure  that  it  possessed  at  hirth. 

Outside  the  hot  spot  (denoted  by  R  in  Figure  1),  i.e.  for  £  =  0(1)  and  beyond, 

H 

the  solution  remains  essentially  frozen  at  its  state  at  the  end  of  the  ignition  stage. 

This  is  because  these  diffusion-dominated  regions  are  unable  to  respond  to  exponentially 
rapid  variations.  Thus,  in  the  inert  region  R  ,  the  solution  is 
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(5.  A) 


0=6  (x,t)  =  e  (x,t  +  £2T  )  *  est 
I  I  C  09 

while  in  R^  and  R^  the  solution  is,  respectively, 

6  =  0,.(/£  X,  t  +  £ 2 T  )  +  Ell'  (X,0)  +  o(£) 
X  c  ®  n 

and 


(S.7) 


9  -  0  (e£,t  +  g2t  )  -  7  £  In  e  +  e  fry  (£)  -  In  A  1  ♦  of  e  )  .  C^.®) 

I  c  »  2  in  n 

The  last  two  expressions  follow  from  (4.1)  and  (4.13),  respectively,  in  conjunction  with  the 
relation  0  =  9j  +  ^. 

The  hot-spot  solution  (5.S)  does  not  match  with  the  P  -  solution  (5.R).  As  pointed 

out  in  (1],  this  is  due  to  the  loaarithmic  behavior  of  t)  (O  for  £  *  n  (cf.  (4.7)), 

which  causes  (5.P)  to  break  down  before  the  edcje  of  R  is  reached.  The  situation  is 

H 

remedied  by  recognizing  that  an  intermediate  region  R^  has  now  evolved  (see  Fiqure  1), 

separating  the  rapidly  shrinking  P  from  P  ,  and  character ized  by  the  variables  o 

H  1 

and  z,  where 

z  =  -£  ln(/£  Aft£2).  0  <  z  <  o  .  (5*9) 

By  rewriting  (2.9)  in  terms  of  the  new  variables,  it  can  he  shown  that  the 
intermediate  expansion  proceeds  as 


c 


(0  -2) 
c 


1  2  -2 

+  c0  (0  -z)  '[an 
c  c  0 


(9 


In 


-z)2{0  +P-e2(8  -z)'1} 
c  c  c  , 


P0  ■ 


(5.  10) 


which  matches  with  the  outer  solution  (5.P)  as  z  ♦  0  and  with  the  hot  spot  solution  (^.S) 
as  z  ♦  o.  This  completes  the  description  of  the  explosion  stage. 


I 
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6.  DETACHMENT  STAGE 


The  hot  spot  reaches  maturity  at  o  =  B0  /(0  +3)*  Then,  as  shown  in  Figure  1, 

c  c 

approaches  its  maxiimsn  value  6^  +  B  and  the  mass  fraction  of  the  reactant,  Y,  * 

zero  inside  the  hot  spot*  At  the  same  time,  a  logarithmic  singularity  develops  in  the 
second  term  of  the  hot-spot  solution  (5.5),  signalling  a  breakdown,  while  the  expansions 
(5*6-5. 8)  and  (5*10),  pertaining  to  regions  outside  the  hot  spot,  still  hold.  Further 

evolution  of  the  hot  spot  is  now  governed  by  the  independent  variables  x  and  s,  where 

2  2 

x  m  6y  and  t=t  +  e  T  +  6  s  , 

A  c  oo 

and  the  exponentially  small  parameter  6  is  defined  by 

62  =  (B/A„  )£1/2exp[-Be  /{2e  (8  +  b)}]  . 

o  c  c 

Since  the  above  scaling  implies 

a  *  B6  /(8  +  B)  -  e  ln( -Bs/e )  , 

c  c 

we  expect  the  detachment  stage  to  match  with  the  explosion  stage  as  s  ♦  In  the  new 

•  •  ■  -  -  —  •  « 

variables,  (2.9)  reads 


8=6  +  (8  +  B  -  8  +  est)exp[ (82/e ) { (8  +  F)-1  -  8-1}]  . 

s  XX  c  c  c 

which  suggests  the  following  expansion: 

8  =  8c  +  B  +  Ew1 (x,8)  +  ...  , 

where 


(6  *  1 ) 


*  2  2 
e  =  <8  ♦  b)  e/8  . 

c  c 

The  problem  for  w  ^  is  found  to  be 

W1  ■  wi  -  w,e  ,  X  >  0.  s  >  -  -  , 
s  XX 


(8.2a) 


w1  (0,9)  =  0  , 
X 


(6.2b) 


-In  x2  “  ln( In  x2) 


(6.  2c) 


w,  ♦  -ln(-s)  -  ln(ln(-s)}  +  90(x/^-»)  -*s  s  +  -  ®,  x/’1'-8  fixed  .  (6. 2d) 


The  argument  behind  the  left  boundary  condition  (6.2b)is  the  same  as  that  behind  (5.4). 
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analogous  to  those  for  the  transition-stage  problem  encountered  in  [1],  and  were  solved  hv 
a  similar  numerical  procedure.  The  solution,  displayed  in  Figure  4,  exhibits  an  expanding 
region  adjacent  to  the  surface  in  which  w^  approaches  zero  as  s  increases.  At  the  same 
time  the  region  of  significant  spatial  structure  is  pushed  farther  and  farther  into  the 
interior.  Thus,  as  the  last  remnants  of  the  reactant  are  consumed,  the  hot  spot  evolves 
into  a  travelling  front,  or  flame,  or  deflagration  wave. 
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Figure  4.  Numerical  solution  of ’  the  detacYiment  stage. 


7 


PPOPA CATION’  STAGE 


In  this  section  the  speed  and  structure  of  the  pronaaa  1  na  wave  are  determined.  The 
structure  is  found  to  be  quasi-steady  and  self-similar  (see  ^inures  1  and  * ) .  It  consists 
of  a  thin  reaction  zone  surrounded  by  a  sliahtly  thicker  envelope  which  separates  the  coll 
material  ahead  of  the  front  fron  the  burnt  material  behind  ( cf ,  '1’  1.  The  portion  of 
envelope  ahead  of  the  reaction  zone  is  a  preheat  reaion  where  the  material  is  brouah*-  up  vo 
the  reacting  temperature. 

Figure  1  shows  that  the  termperature  of  the  cold  material  into  which  the  front 
advances  falls  steadily  awav  from  the  surface.  This  leads  to  a  oradual  retardation  an’ 
thickening  of  the  flame,  as  we  shall  see.  Let  x0,  taken  to  be  0(1),  be  the  tine- 
varying  location  of  the  flame.  We  shift  to  a  reference  frame  in  which  the  flame  is  at  rest 
by  setting 

X  =  x0  +  (S  ,/c  .  0.1) 

where  c  is  the  spatial  coordinate  in  the  envelope.  The  small  parameters  >5^ 
(exponentially  small  in  e)  and  (proportional  to  c),  themselves  dependent  upon 

Xg,  are  to  be  assigned.  The  scalinq  (7.1)  clearly  implies  that  6^/c^  is  the  thickness 
of  the  envelope.  We  shall  later  find  that  6^  is  the  thickness  of  the  react  ion  zone 
within  the  envelope.  The  appropriate  time  scale  r  for  the  propaaation  staae  is  defined 
via  the  expression 

t  =  t  +  e2T  *  (6  /t  )  r ,  r  >  0  ,  H.2) 

C  °°  11 

indicating  that  the  solution  continues  to  envolve  with  exponential  rapidity.  The  wave 
speed  now  emerges  as 


% 

dt 


% 

dr 


t  U(r) 

S 


(7.1) 


say,  where  U(r)  is  0(1).  The  inert  temperature  at  Xq  will  play  a  critical  role? 
therefore  we  set 

9°  =  6  (X  ,t  +  E2T  )  (7.4) 

I  I  n  c  00 

and  observe,  from  (2.7)  and  (7.2),  that  for  r  =  0(1), 
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(7.5) 


8  (x  ,t>  =  6j  +  est  • 

we  can  now  define  61  and  E  ^  as  follows: 

«2  =  (B/A„)e1/2exp!(62/£){(6°  +  B)'1  -  0'1}]  (7.6a) 

in  cl  c 

and 


e ,  =  (6°  +  B)2/e  =  (e"  +  B)2£/02  . 

1  I  I  C 

Thus,  both  6^  and  e ^  depend  upon  xn  through  0^.  Since  6^ 
(to  urity)  as  Xg  increases  (cf.  (2.7)  and  (7.4);  also  see  Figure 
decreases  with  Xq.  However,  the  ratio  6 ^/e  increases,  so 
advances  it  thickens  and,  from  (7.3),  its  speed  falls. 

In  the  C,r  variables,  (2.9)  transforms  into 

61(e1  -  rU(30"/3xn)]0r  =  8  +  u8^  + 

e~2(0i  +  B  -  0)exp[{(8°i  +  B)2/e}{(0°  +  B)_1  - 


(7.*b) 

decreases  monoton ical lv 
1),  6^  increases  while 

that  as  the  flame 


0-1}]  ,  <7-7) 
r  >  0  . 


The  subsequent  analysis  follows  along  the  lines  of  that  in  section  6  of  ft],  so  that  a 

brief  account  should  suffice.  Separate  solutions  need  to  be  constructed  for  the  reaction 

zone,  the  envelope  and  the  region  outside  the  front.  In  this  outer  region  (see  Figure  1), 

(  0  (x,t)  +  B  =  0  (x,t  +  e2T  )  +  B  +  est,  x  <  0  (burnt  reqion)  ,  (7. “a) 

x  I  c  00  0 


0  = 


l  0j(x,t)  =  BjfXft^  +  e2^)  +  est,  Xq  >  0  (cold  region)  .  (7.Rb) 

The  burnt-region  solution  (7.fla)  is  an  exact  solution  of  (7.8),  satisfying  the  surface 
boundary  condition  in  (2.2a),  while  the  cold-region  solution  (7.8b)  is  just  the  frozen 
solution  (5.6).  In  the  envelope  the  portion  £  <  0  is  a  part  of  the  burnt  region,  i.e. 

0  =  0°  +  B  +  est,  5  <  0  .  (7 . da) 

For  5  >  0,  chemistry  is  still  unimportant  since  0  <  0°  +  B.  The  leading-order  balance 
in  (7.7)  is  therefore  diffusive-convective,  yielding 

0  =  A  ^ ( s)  +  B^t  s) exp  [-U( s)£]  +  ...  .  (7.9b) 

Continuity  of  (7.9a,b)  across  £  ”  0  gives 
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*1  +  pi  -  6I  +  P 


while  matching  between  the  preheat  solution  (7.9b)  and  the  cold-region  solution  (7,Pb) 
leads  to 

A,  =  e°X  *  <7' 

Then,  (7.10)  determines  i.e. 


The  remaining  unknown  U  is  determined  by  analyzing  the  structure  of  the  reaction  zone. 
This  is  accomplished  by  the  reaction-zone  scaling 

C  =  e,p 

(which  confirms  that  6  is  the  thickness  of  the  reaction  zone)  and  the  reaction-zone 


expansion 


0  =  +  B  +  t^vlp.r)  +  ...  , 

which  reduce  (7.7),  to  leading  order,  to 

v  "  v  exP  v=0,  -*»  <  p  <  “  . 
PP 

Matching  with  (7.9a,b)  provides  the  boundary  conditions 

v+0  as  p  •*  -  «,  v  ♦  -  B1Up  as  p 
which,  when  applied  to  the  first  integral  of  (7.13),  yield 


whence  (7.12)  leads  to 


8^-/2 


U(s)  =  /2/B  . 


substituting  for  0,  e1  and  { 1  in  (7.3),  the  speed  of  the  flame  is  seen  to  be 


dx  (0°  +  B)2  /2A  c  f@2 

_ 0  I _  / 0  c 


:r exp  if  -  To — : +  0 


U. 


Tbe  ultimate,  steady  flame  speed  will  correspond  to  0“  »  1  in  the  above  expression;  the 
right  ineouality  of  (2.5)  guarantees  that  this  speed  will  not  be  exponentially  small. 

•Hie  above  analysis  is  predicated  upon  x0  being  0(1),  which  leads  to  the 
temperature  immediately  ahead  of  the  flame  being  essentially  due  to  inert  heating.  Such  is 
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not  the  case  when  the  flame  is  verv  clos^  to  the  surface.  let  xn  he  exponential  lv  s^all, 
i.e.  let 


xn  =  U3/2/V*  "  ' 

Then#  the  front  edoe  of  the  flame  lies  in  the  intermediate  reaion  disrussed  in  faction  S, 

and  there  the  temperature  differs  from  the  inert  solution  by  an  0(1)  amount  ( cf . 

(5.10))*  It  can  be  shown  that  (7.14)  must  then  be  replaced  by 

U(s)  =  /2/ie  +  B  -  6 2 ( 0  -  z„>’  . 

c  c  c  n 

As  z.  0,  (7.14)  is  recovered.  As  z  ♦  B9  /9  ♦  R)  ,  t?  becomes  singular,  snaoest  i  nu  a 

0  0  c  c 

breakdown  of  the  analysis.  This  is  to  be  epxected ,  because  the  flame#  at  that  time#  is 
still  in  the  detachment  staqe  discussed  earlier. 
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8.  CONCLUDING  REMARKS 


An  asymptotic  description  of  the  evolution  of  a  deflaqration  wave  in  an  initial  lv 
material  under  the  stimulus  of  a  constant  heat  flux  has  been  given.  The  analysis  revelas 
that  the  situation  is  quite  similar  to  that  for  self- induced  combusion,  discussed  earlier 
[1],  with  the  exception  that  the  thermal  explosion  process  responsible  for  the  birth  of 
hot  spot  is  now  confined  to  a  thin  surface  layer.  Furthermore,  it  is  found  that  once 
ignition  has  occurred,  events  evolve  so  rapidly  that  the  system  hecomes  oblivious  to  t^p 
presence  of  the  external  heat  supply. 

To  be  sure,  the  results  are  based  upon  the  assumption  of  the  Lewis  number  beinq 
unity.  However,  the  case  L  =  «>,  corresponding  to  solid  combustibles,  was  also 
analyzed.  The  full  system  (2. 1-2.3)  must  then  be  dealt  with  head  on,  since  (2.8)  no  longer 
obtains.  The  reactant  equation  (2.1b)  does  simplify  due  to  the  absence  of  the  diffusion 
term,  and  the  boundary  conditions  (2.2b)  have  to  be  abandoned.  Otherwise,  the  analysis 
proceeds  in  an  analogous  way.  Up  to  the  end  of  the  explosion  period,  such  calculations  ror 
the  case  of  self-induced  combustion  appear  in  the  forthcoming  paper  by  Kassoy  and  Poland 
[3] . 

The  results  differ  from  the  L  »  1  case  only  in  one  significant  respect.  As  lono  as 
combustion  is  confined  to  a  stationary  sone  adjacent  *0  the  surface  ( i  .e  prior  to  the 
detachment  of  the  hot  spot),  any  reactant  consumed  there  is  not  replenished  by  diffusion 
from  outside.  Therefore,  for  a  given  temperature  rise,  the  reactant  concentration  falls  to 
a  lower  value  within  the  zone  than  it  would  if  diffusion  of  reactant  were  permitted.  This 
effect  appears  as  early  as  the  transient  stage.  Thus,  Linan  and  Williams  [11)  found  that 
an  0(e)  change  in  4  corresponds  to  a  (much  larger)  Ol^e )  change  in  Y.  The  same 
holds  true  during  the  ignition  period.  By  the  time  the  explosion  stage  ends,  it  is  seen 
(3]  that  the  maximum  temperature  reached  in  the  stationary  hot  spot  is  lower,  by  as  much  as 
an  0(1)  amound,  than  the  corresponding  value  for  the  case  L  -  1.  The  propagatinq  flame 
does  not  suffer  from  such  a  fate  because  _it  obtains  the  necessary  reactant  supply  throuqh 
convection. 
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